Trapped Fermions across a Feshbach resonance with population imbalance 
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We investigate the phase separation of resonantly interacting fermions in a trap with imbalanced 
spin populations, both at zero and at finite temperatures. We directly minimize the thermodynam- 
ical potential under the local density approximation instead of using the gap equation, as the latter 
may give unstable solutions. On the BEC side of the resonance, one may cross three different phases 
from the trap center to the edge; while on the BCS side or at resonance, typically only two phases 
show up. We compare our results with the recent experiment, and the agreement is remarkable. 
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The recent experimental realizations of condensation 
and superfluidity in Fermi gases near a Feshbach reso- 
nance have provided a new platform for controlled study 
of many-body physics in the strongly correlated region 
0. A significant advance in this direction is illustrated 
by two very recent experiments 0, j which study the 
fermi superfluidity with controlled population imbalance 
of different spin components. With increasing population 
imbalance, a phase transition from the superfluid to the 
normal state has been observed |2| , and phase separation 
of the Fermi gas in the trap has been identified [2, y| • 

The Cooper pairing in the Fermi gas, which is essential 
for the Fermi condensate and the superfluidity, requires 
an equal number of atoms from both spin components. If 
the atomic gas has imbalanced population in its two spin 
states, inevitably some of the atoms will be unpaired, 
which triggers the competition between the Cooper pair- 
ing and the population imbalance. Such a competition 
may lead to various new phenomena Jj J a J P' |8( , which 
have raised strong interest recently [fj L3] . 

In this paper, we present our study on the trapped 
fermions across a Feshbach resonance with population 
imbalance. Compared with other recent theoretical 
works 0, Q, 0, we demonstrate the following results: 
(i) We consider a wide resonance with parameters cor- 
responding to the recent 6 Li experiment, and calculate 
the pairing gap and the density distribution, which can 
be directly compared with the experiment, (ii) We take 
into account the trapping potential, and show that the 
trap favors the phase separation, which is critical for the 
understanding of the recent experiments 0, Q Some 
of the co-existing or unstable phases in a homogeneous 
gas @ necessarily become phase separated even if a very 
weak trap is turned on. (iii) We investigate the effect 
of the population imbalance both at zero and at finite 
temperatures. At zero temperature, the BCS superfluid 
state requires an equal number of atoms from both spin 
components; but at finite temperature T it allows popu- 
lation imbalance carried by the quasiparticle excitations. 
This key difference significantly relaxes the competition 
between the population imbalance and the Cooper pair- 
ing at finite T, which helps to explain the robustness of 



the superfluid state in the recent experiment 0. (iv) We 
directly minimize the thermodynamical potential to find 
out the stable configuration of the system, and caution 
the use of the gap equation in the case of population im- 
balance. In the region with competition between different 
phases, the gap equation often gives incorrect results or 
unstable solutions. 

To describe the fermions across a Feshbach resonance, 
we take the following two-channel Hamiltonian [9I IT0II1H : 
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where ek = k 2 /(2m) (m is the atom mass and h = 1), 
fia is the chemical potential for the spin-er component 
(a =t,| labels the two spin states), V is the quanti- 
zation volume, a and & q are the creation operators 
for the fermionic atoms (open channel) and the bosonic 
molecules (closed channel) , respectively. The bare atom- 
molecule coupling rate a, the bare background scatter- 
ing rate U, and the bare detuning 7 are connected with 
the physical ones a p , U„, j p through the standard renor- 
malization relations 11], and the values of a p , U p ,j p are 
determined from the scattering parameters (see, e.g., the 
explicit expressions in Ref. We take the local den- 

sity approximation so that /if = /i r + h, fi± = /i r — h, 
fi T = fi—V(r), where V (r) is the external trap potential 
(slowly varying in r), and [i,h are determined from the 
total atom number N — jVj + N± and the population im- 
balance /3 = |iVj — jVjJ /N through the number equations 
below. 

For simplicity, we first take a mean-field 
approach by assuming (6 q ) = (&o)<5q0 = 
— a/yV X) k ( a -k.i a k,|) / (7 ~ 2/x r ) (the last equality 
comes from the Heisenberg equation for the operator fe q ). 
At this level, we neglect the pair/molecule fluctuation, 
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an approximation that is valid near zero temperature 
or in the BCS limit. Later we will discuss how to 
incorporate the pair/molecule fluctuation directly into 
the final equations. We also neglect here the possibility 
of a non-zero-momentum pairing (the FFLO state [4|, 
with the pair momentum q ^0). This is motivated by 
the fact that the FFLO state is stable only within a 
narrow parameter window deep in the BCS region ||| 
and is absent in the recent e Li experiments 0, 0] • 

With the mean-field approximation, one can 
then find out the thermodynamical potential 
il = — Tln[tr(e~ ff / T )] at temperature T (the Boltzmann 
constant is taken to be 1). It is given by the expression 

n = -A 2 V/U T - T]T{ln[I + exp(- \E n \ /T)] 

k 

+ ln[l + exp(-£ kT /T)] - [e k - H - 9(E kl )E kl ]/T}, (2) 

where the parameter Ut = U — a 2 / (7 — 2/i r ), the quasi- 
particle excitation energy E k -\^ = "J ( e k — Hr) 2 + |A| 2 =p 
h, and the gap A = a{b )/VV + U/V Ek( a -k,l a k,T> = 
z (bo) I VV (z = a — U (7 — 2/i r ) / a) . The ^-function is de- 
fined as 9 (x) = 1 for x > and 9 (x) = otherwise. Note 
that different from the equal-population case, the excita- 
tion energies E ka are different for the a — t, | branches. 
Without loss of generality, we take h < so that E k ^ > 
always; while for a certain range of h, the sign of E k ± de- 
pends on the momentum k . When E k i < 0, the atoms 
at that momentum are actually not paired (and thus have 
no contribution to A) as pairing is energetically unfavor- 
able (that is why we have \E k i\ and 9(E k i) in Eq. (2)). 
This case corresponds to the so-called breached pair state 
[l3|. From the variational condition dfl/d\A\ 2 — 0, we 
get the gap equation 
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where the parameter 1/Ut p — 

V [Up ~ ay ( lp - 2 Mr )] = l/U T + (1/V) Ek [1/ (2*)] 
(the latter equality comes from the renormalization 
relation between j,a,U and j p ,a P) U p ^3). From the 
relation dQ/d/-i a = —n ra V, where n ra is the density of 
the spin-cr component, we get the number equations 
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[ulf(E Ka ) + v 2 f(-E K ^)} 
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where the parameters uf. = (E k + (e k — fi r ))/2Ek, v£ = 



(E k -(eu-fr))/2E k , E k = y/(e k - n r ) 2 + |A| Z , thefcrmi 
distribution f(E) = 1/ (l + e E / T ), and we take - T=4 
and vice versa. The last term |A| / z 2 m n rcr comes 
from the contribution of the molecule (closed-channel) 
population, which is small near a wide resonance. The 
atom densities n r | and n r | are connected with the total 



atom number and the population imbalance through iV = 
/ d 3 r (n r f + n r |), and /3 = J d 3 r5n r /N (6n r = n r |— n r |). 

The mean-field approach above neglects the 
pair/molecule fluctuation. With this fluctuation taken 
into account, there will be a non-condensed fraction 
of the pairs, which also contributes to the gap for the 
fermionic quasiparticles S° the gap now is replaced 
by |A| 2 = |A S | 2 + |A p J^ where |A S | 2 = z 2 |(6 )| 2 /V 
and |A P0 | = z 2 m nc /V with |(&o)| and m nc repre- 
senting the condensed and the non-condensed molecule 
numbers, respectively. With the interpretation of |A| as 
the total gap, the gap and the number equations above 
remain valid with the pair fluctuation (note that the 
contribution to n ra from the non-condensed molecules 
m nc is automatically taken into account by the last term 
in Eqs. (4)). If one wants to break up A to find out 
the superfluid order parameter A s and the pseudogap 
A pg , one needs to know the dispersion relation for the 
pair/molecule excitations. Here, to compare with the 
experiments 0i , we only calculate the total gap A and 
the atom density distributions n rcr . For that purpose, 
the above gap and number equations suffice, and we do 
not need to specify the pair dispersion relation. Note 
that the total gap A is the quantity directly measurable 
through the radio- frequency spectroscopy |l5| . 

The gap and the number equations above are in prin- 
ciple sufficient to determine the distribution of |A| and 
n rcr . However, different from the equal-population case, 
the solution of this set of equations turns out to be subtle, 
as they often give unstable phases or incorrect results. To 
understand this, we examine the behavior of the thermo- 
dynamical potential fl under the variations of some sys- 
tem parameters. Note that the gap equation should give 
an optimal value of |A| which minimizes Q. However, 
with a population imbalance, the potential f2 has a dou- 
ble well structure in many cases. In Fig. I, we show f2 as 
a function of |A| as we change the chemical potential /i r 
(from the trap center to the edge), the potential differ- 
ence h (with varying population imbalance), the system 
temperature T, and the field detuning "f p . One can see 
that the double- well structure shows up, in particular on 
the BCS side or in the low-temperature, large-imbalance 
cases. Typically, one minimum of the double-wells corre- 
sponds to |A| =0 (the normal state) and the other has 
a nonzero |A|. When has a double- well structure, the 
gap equation can give an incorrect result in several differ- 
ent ways: first, it may pick up a solution corresponding 
to the maximum of Q (the gap equation is satisfied at 
that point), which is obviously not stable; second, it can 
choose the shallow well, which gives a metastable but not 
the optimal configuration; and finally, the gap equation 
(3) is usually not satisfied with |A| = even if the latter 
is the minimum of (as dfl/d \ A\" ^ at that point). In 
this case, the numerical program gives an incorrect result 
in an uncontrollable way. 
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FIG. 1: The thermodynamical potential fl as a function 
of the total gap |A| with varying (a) h (the chemical po- 
tential difference), (b) 7 P (the field detuning), (c) \i (the 
chemical potential), and (d) T (the temperature). We have 
used Tf = Ef = k F /2m as the unit of energy/tempeature, 
where &f = (37r 2 no) 1,/3 is a convenient inverse length scale 
corresponding to a density no = 3 x 10 13 cm -3 , typical for 
the 6 Li experiments. The parameters a p , U p take the stan- 
dard values for the 6 Li atoms with a vy friQ — — 0.6Ef, and 
U p no ~ —82Ef- The other parameters for Figs, (a-d) are 
given by (a)u = 0.5E F , 7 = 0, T = 0.01T F ; (b) fi = 0.5E F , 
h = -0.35.Ef, T = 0.01T F ; (c) h = -0.35.Ef, 7 = 0, 
T = O.OITf; and (d) fi = 0.5E F , h = -0.35.Ef, 7 = 0. 



To overcome the above problems, we always check the 
stability of the solution by finding out all the minima of 
the thermodynamical potential f2. Note that although 
the metastable state at the bottom of a shallow well of f2 
does not give the ground state configuration, it has a fi- 
nite relaxation time and may appear in real experiments 
under certain circumstances (similar to the superheating 
or supercooling phenomena in a classical phase transi- 
tion). So, in that sense, the state of the system may 
not be unique. To avoid this complexity, in the following 
we assume the system to be in its equilibrium configura- 
tion by choosing the solution corresponding to the global 
minimum of the thermodynamical potential Q. 

Through minimization of the thermodynamical poten- 
tial, we have calculated the gap and the density distribu- 
tions for trapped fermions with different magnetic field 
detunings, at both zero and finite temperatures. Fig. 2 
shows some typical results. To summarize, we find that 
the gas is separated into several different phases from the 
trap center to the edge, with the number of phases and 
their boundary sensitive to the detuning j p , the popula- 
tion imbalance /3, and the temperature T. At zero tem- 
perature and on the BEC side of the resonance with a 
small population imbalance, there are typically three dif- 
ferent phases: the superfluid phase occupies the center of 
the trap, where all the particles are paired with no pop- 
ulation imbalance (the SF phase). Further out from the 



FIG. 2: The gap (|A|) and the differential density (|A| and 
5n r ) distributions for the trapped fermions, both at zero and 
at finite temperature. The x-axis is the distance from the 
trap center, in the unit of iyr. The solid curves represent 
|A| (the left j/-axis) and the dashed ones are for Sn r ( the 
right y-axis). The other parameters are given by (a)(c)/x = 
-111.48.Ef, h = -111.52.Ef; and (b)(d) u = 0.5E F , h = 
—O.IEf- The corresponding population imbalance ratio 
and the total atom number N are (a) 1%, 3.1 x 10 5 ; (b) 
6%, 3.9 x 10 7 ; (c) 69%, 7.9 x 10 5 ; and (d) 17%, 4.0 x 10 7 , 
respectively. 



center, when the pairing gap becomes smaller than the 
chemical potential difference \h\, there is a second phase, 
in which the particles are paired except for a finite re- 
gion of the momentum shell (the breached pair phase [|J , 
or in short, the BP phase). This BP phase is charac- 
terized by a nonzero |A|, but with excessive fermions in 
the breached momentum shell, which gives a finite pop- 
ulation imbalance. In Fig. 2(a), the middle region with 
|A| and Sn r both nonzero corresponds to such a phase. 
Towards the edge of the trap, where the pairing gap has 
already vanished, only the normal Fermi sea is left, with 
different Fermi surfaces for different spin components. As 
the population imbalance increases, the SF phase at the 
center of the trap depletes and yields to the surrounding 
BP phase and the normal Fermi sea. At the resonance 
or further out to the BCS side, this picture is different in 
that there are now only two stable phases, the SF phase 
and the normal phase, as shown in Fig. 2(b) (the BP 
phase in the middle loses its stability). As the popula- 
tion imbalance increases and reaches a critical value, the 
system undergoes a SF-normal phase transition, and the 
trap is left with only the normal Fermi gas. 

At finite temperature, the main difference is that there 
are fermionic excitations in the SF phase, which carry 
population imbalance of the spin components. As a re- 
sult, it becomes easier to satisfy the population imbal- 
ance constraint, as with a fixed imbalance ratio /?, the 
corresponding chemical potential difference \h\ becomes 
significantly smaller. This helps to stabilize the SF phase 
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come the instability problem associated with the solution 
of the gap equation. We establish a general phase sepa- 
ration picture for the trapped fermions across the whole 
resonance region, at both zero and finite temperatures. 
We have also compared our calculation with a recent ex- 
periment, and recovered some of the main experimental 
findings. 

Note added: Upon completion of this work, two 
preprints [P. Pieri, G.C. Strinati, |cond-m at/0512354 J. 
Kinnunen, L.M. Jensen, P. Torma, cond-mat/05f 2556 
appeared on arxiv, where a similar problem are investi- 
gated at zero temperature with different theoretical ap- 
proaches. 



FIG. 3: The atom density distributions in a trap, calculated 
with parameters corresponding to the recent MIT experiment 
with T ~ 0.1E F , and 7 ~ 2100£f {B - B ~ 56.3G). 
The solid, the dashed, and the dash-dot curves correspond to 
the imbalance ratio /3 ~ 46% (the corresponding fi = 0.5Ef, 
h = -0.12.E F ), (3 ~ 86% (the corresponding /i = 0A3E F , 
h = -0.23.Ef), and /3 = (fj, = 0.5E F , h = 0) respectively, 
with the total atom number N ~ 2.7 x 10 7 in all the cases. 
Figs. (c,d) show the column integrated density distributions 
corresponding to Figs. (a,b) . 

in the case of imbalanced population. Because of this fea- 
ture, we cannot use Sn r to distinguish the SF and the BP 
phases any more, so the boundary between them becomes 
obscure. The population difference Sn r is peaked at the 
point where the gap |A| vanishes, marking the only dis- 
tinguishable phase separation between the paired phase 
(SF or BP) and the normal phase (see Fig. 3c and 3d). At 
higher temperatures, the paired phase in the trap shrinks 
and finally disappears at the pair disassociation temper- 
ature T* (the superfluidity should disappear before that 
with T c < T* 11]). In Fig. 2(c) and (d), we show the 
distribution of |A| and 5n r at a typical temperature of 
0.1-Ef, for the BEC and the resonance-BCS sides, respec- 
tively. 

Finally, to compare with the recent MIT experiment 
, we calculate the density distributions n ra , with all 
the parameters roughly the same as those in the experi- 
ment (the experimental temperature may have some un- 
certainties). The results are plotted in Fig. 3. It shows a 
pretty good agreement with the experimental data (Fig. 
3 of Ref. 0), at least semi-quantitatively. In particular, 
the calculation shows that the distribution 5n r has a dip 
at the trap center with a moderate population imbalance 
((3 = 46%); while the dip disappears when (3 becomes 
large (86%), which agrees exactly with the experimental 
findings. 

In summary, we have studied the effects of trap and 
temperature on fermions across a wide Feshbach reso- 
nance with population imbalance. We propose to directly 
minimize the thermodynamical potential in order to over- 
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